Draft version March 1, 2013 

Preprint typeset using X^T^^X style emulatcapj v. 5/2/11 



(N 

o 

(N 



(N 
O 



CONSTRAINING TYPE la SUPERNOVA MODELS: SN 201 Ife AS A TEST CASE 
F. K. Ropke/'^ M. Kromer,^ I. R. Seitenzahl/^ R. Pakmor,' S. A. Sim,^ S. Taubenberger,^ 



F. ClARALDI-SCHOOLMANN,- W. HiLLEBRANDT," G. AlDERING, 



Fakhouri,'"' M. Fink/ 
M. Kerschhaggl,' 

.9 



^ D. 

M. 



Fouchez, E. Gangler, J 



P. Antilogus,'' C. Baltay,' S. Benitez-Herrera/ 
M. Childress,^'^ N. Chotard,"' Y. Copin,"' H. K. 



Guy," S. Hachinger,^ E. Y. Hsiao,'' C. Juncheng," 



S. Perlmutter,' D. Rabinowitz 



KowALSKi, p. Nugent, K. Paech, R. Pain, E. Pecontal 

C. Saunders, ■'''''' G. Smadja,^" 

13 



10 



M. RiGAULT, 

R 



K. RUNGE,' 



R. Pereira 
N. Suzuki,'^ C. Tag 



13,14 



C. Thomas," A. Tilquin," G. Wu"-^'^ 

Draft version March 1, 2013 

ABSTRACT 

The nearby supernova SN 2011fe can be observed in unprecedented detail. Therefore, it is an im- 
portant test case for Type la supernova (SN la) models, which may bring us closer to understanding 
the physical nature of these objects. Here, we explore how available and expected future observa- 
tions of SN 201 Ife can be used to constrain SN la explosion scenarios. We base our discussion on 
three-dimensional simulations of a delayed detonation in a Chandrasekhar-mass white dwarf and of a 
violent merger of two white dwarfs — realizations of explosion models appropriate for two of the most 
widely-discussed progenitor channels that may give rise to SNe la. Although both models have their 
shortcomings in reproducing details of the early and near- maximum spectra of SN 201 Ife obtained 
by the Nearby Supernova Factory (SNfactory), the overall match with the observations is reasonable. 
The level of agreement is slightly better for the merger, in particular around maximum, but a clear 
preference for one model over the other is still not justified. Observations at late epochs, however, 
hold promise for discriminating the explosion scenarios in a straightforward way, as a nucleosynthesis 
effect leads to differences in the ^^Co production. SN 2011fe is close enough to be followed sufficiently 
long to study this effect. 

Subject headings: Supcrnovae: general — supernovae: individual (SN 2011fe) — hydrodynamics — 
nuclear reactions, nucleosynthesis, abundances 
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Perhaps the most fundamental problem hindering a 
better understanding of SN la explosions is the unclear 
nature of the progenitor system. One way of addressing 
this problem is to carry out numerical simulations for 
different scenarios that involve thermonuclear explosions 
of white dwarfs (WDs) and to compare the results with 
observations. Obviously, detailed observational data are 
a prerequisite for this approach. At the same time, the 
comparison should be based on models that avoid free 
parameters in the description of the explosion mecha- 
nism, as far as possible. Only then the predictive power 
of theoretical models will be sufficient to discriminate be- 
tween explosion models and to draw conclusions about 
progenitor systems. 

In addition to the possibility of directl y con straining 
the progenitor system from archival d ata (iLi et al. 20li 
Liu et al.ll201lD or early observations ([Brown et al.ll2011 



Nugent et al.ll2011bl:lBloom et al.|[2012[) . the recentlv dis 
covered nearby SN la 201 Ife offers a unique opportunity 
for a comparison with explosion models. Of particular 
value are the possibility to follow this close object photo- 
metrically to extremely late epochs and the exact knowl- 
edge of the explosion time. 

SN 201 Ifc was first detected by the Palo mar Transient 
Factor y on 2011 August 24.167 in MlOl (jNugent et al.l 
I2011a[) . A preliminary analysis of our data indicates 
that it reached an apparent _B-band peak magnitude 
of 9.9 on September 11. Combined with the de- 
rived e xplosion date of 2011 August 23.7 (jNugent et al.l 
l2011bD . the S-band rise time of SN 201 Ife is ~18 3d— 
a typical value for normal SNe la (jConlev et al.l 120061 : 
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iHavden et al.ll2010D. Assuming a distance to MlOl of 



6.4 Mdc (IShappee fc Stanea[20Tll) . SN 201 Ife is a normal 
SN la with Ms^max = —19.13, having produced ^--^0.6 Mq 
of 5^Ni (|Stritzinger et al.ll2006l ). The identification of 
SN 20 life as a prototypical SN la is also corroborated 
by the observed spectra (as shown below). 

With the development of three-dimensional simula- 
tions of thermonuclear explosions in carbon-oxygen WDs 
and of the subsequent radiative transfer (RT) leading 
to the formation of the observables, a new generation 
of models is currently becoming available. These have 
the advantage that the explosion physics is represented 
in a far less parameterized manner than in previous 
one-dimensional models. Due to their improved predic- 
tive power, a comparison with observational data would 
in principle allow us to constrain the explosion sce- 
nario of SNe la. However, no currently available multi- 
dimensional model reaches the level of agreement that 
can be obtained fitting one-dimensional semi-empirical 
models to data. This challenges the interpretation of the 
comparison between the new models and SN la data. 

Here, we address the question of whether SN 201 Ife can 
be explained by models of an exploding Chandrasekhar- 
mass WD (realized as a delayed detonation) or a violent 
merger of two WDs. These scenario s can lead to observ- 
ables that resemb l e normal SNe la (iMazzali et al.l 120071 : 
iKasen et al.1 120091: iPakmor et all 12010), but they differ 
fundamentally in the explosion mechanism, the mass and 
the structure of the ejecta. A discrimination between 
them based on comparison with observations would help 
to shed light on the open question of the progenitor sys- 
tem. An explosion of the WD near the Chandrasekhar 
mass is usually attributed to the single-degenerate pro- 
genitor channel in which a carbon-oxygen WD accretes 
matter from a non-degenerate companion; however, the 
formation of a Chandrasekhar-mass object in a merger 
of two WDs cannot be excluded. Our second scenario 
results from the merger of two WDs with similar and 
rather high masses adding up to a total of 2 AIq. Both 
models are set up to produce ^0.6 Mq of ^^Ni but apart 
from that they follow generic assumptions and are not 
tuned to fit the data of SN 201 Ife. 

2. EXPLOSION MODELS 

The most promising way of producing observables 
in reasonable agreement with observations of normal 
SNe la from an explosion of a Chandr asekhar-mass WD 
is the delayed detonation mechanism (jKhokhlov| |1991|). 
We model this scenario u sing the techniques describe d 
by iReinecke et all (Il99 9l): lR6pke fc Hillebrandtl (l2005l) : 
ISchmidt et al.l (|2006f ) and lRopke fc Niemeveii (|2007[ ). 

An isothermal (T = 5 x 10^ K) WD composed of car- 
bon and oxygen in equal parts by mass was set up in 
hydrostatic equilibrium with a central density of 2.9 x 
10^ g cm~'^ and an electron fraction of Ye = 0.498864, 
corresponding to solar metallicity. The model was dis- 
cretized on a three-dimensional Cartesian moving grid 
(jRopkel [20051 ) with 512"^ cells consisting of two nested 
parts. To reach the intended ^^Ni production, the initial 
deflagration was ignited in 100 sparks placed randomly 
in a Gaussian distribution within a radius of 150 km from 
the WD's center on the inner grid, which had a resolution 
of 1.92 X 10^ cm. After an initial deflagration phase simi- 
lar to that described bv iRopke et al., (,2007. ) . a detonation 



was triggered at every location on the flame where the 
fuel density was in the range of 6 to 7 x 10^ g cm~^ and 
the grid cell contained preferentially fuel material, pro- 
vided that the turbulent velocity fluctuations exceeded 
10* cms~^ at a signiflcant fraction of the flame area and 
persisted for sufficiently long times. This loo sely follows 
the criteria proposed bv lWooslev et al.l ()2009[ ). Since the 
initiation of a detonation proceeds on scales that are not 
resolved in our simulations, the probability of finding 
high turbulent velocities on unresolved scales is extrapo- 
lated applying the procedure of|R6pke (2007). The evo- 
lution was followed to a time of 100 s after ignition, by 
which homologous expansion of the ejecta was reached to 
a good approximation. This model, called NlOO, is part 
of a larger set of delayed-detonation simulations (Seiten- 
zahl et al., in preparation). 

The details of the nucleosynthesis in this explo- 
sion were determined from thermodynamic trajectories 
recorded by 10^ tracer particles distributed in the explod - 
ing WD (jTravaglio et all [200l iSeitenzahl eFaLl [20Tot) . 
The characteristics of the model are summarized in Ta- 
ble ffl 

The second simulation we discuss models the inspi- 
ral, merger and explosion of two WDs with 1.1 Af© and 
0.9 Mq, respectively. Details of the corres ponding sim- 
ulations are given by IPakmor et al.l (|2012D . While the 
inspiral and merger phases were followed with a version 
of the SPH code Gadget (Springel 2001), the subse- 
quent thermonuclear detonation was modeled with tech- 
niques similar to those employed in NlOO. The ques- 
tion of whether a detonation triggers at the interface 
between the two merging stars is controversial. Sim- 
ulations with sufficiently high numbers of SPH parti- 
cles, such as presented here, show the formation of a 
hot spot, and we assume a detonation to trigger at this 
location when the temperature exceeds 2.5 x 10^ K in ma- 
terial of p « 2 X lO^gcm"^, relying on the microscopic 
simulat ions of detonation initiation by ISeitenzahl et al] 
()2009al) . Again, the evolution was followed up to 100s 
and the composition of the ejecta was determined in a 
post-processing step. The results of this simulation are 
given in Table [TJ 

Density and composition of both models in homologous 
expansion are visualized in Figs. [T]and[2] We note that 
the ejecta structure resulting from the WD- WD merger 
differs fundamentally from that of NlOO. This is because 
the explosion of the secondary WD happens shortly af- 
ter that of the primary. Therefore, the outer ejecta ma- 
terial originates from the primary. At the onset of the 
explosion, the primary had a radius below 0.02 Rq mak- 
ing our violent merger scenario consistent with the con- 
straint on the radiu s of the exploding object derived by 
iBloomet all (J20T1 . 

3. COMPARISON WITH SPECTRA OF SN 2011fe 

From the nucleosynthesis tracer particles we con- 
structed detailed abundance distributions of the explo- 
sion ejecta at 100 s and mapped them to 50^ Carte- 
sian grids. These grids were then used to derive syn- 
thetic light curv es and spectra with t he Monte Carlo 
RT code Artis (|Kromer fc SimI l2009l : ISm]l2007l) . To 
this end, we simulated the propagation of 10* photon 
packets from 2 to 120 da ys after ex plosion using the 
cd23_gf-5 atomic dataset of iKromer fc Sim. (,2009l) . which 
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Figure 1. Slices through our delayed-dctonation model NlOO in the x— 2:-plane showing the density (top left) and abundance distribution 
of selected species at 100 s after explosion. 



Table 1 

Model characteristics. 



delayed detonation (NlOO) violent merger 



total ejecta mass [Mq] 

asymptotic kinetic energy of ejecta [10^^ 



s^Ni mass [Mq] 

total iron group [Mq] 

total intermediate mass elements [Mq\ 

carbon mass [Mq] 

oxygen mass [Mq\ 

combined mass of ^^Fe and ^^Co {Mq\ 

combined mass of ^^Ni and ^'^Co [Mq] 



erg] 



1.40 
1.45 



1.95^ 
1.7 ~ 



0.604 

0.839 

0.454 

0.003 

0.101 

1.33 X 

1.88 X 



nSOU — ^ 



£)1000 

late 



: ™-1900d 



- mgood in leptonic light curve [mag] 2.25 

- mgood in leptonic light curve [mag] 3.20 



0.616 

0.697 

0.5 

0.153 

0.492 

3.73 X 10-3 

1.49 X 10-2 



B-band rise time [days] 


16.6 


20.8 


B-band peak luminosity [mag] 


-19.0 


-19.0 


Ami5(-B) [mag] 


1.34 


0.95 



2.65 
3.87 



''0.05 Mq are lost during the explosion simulation because of the finite extent of the grid 
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Figure 2. Same as Fig. [l] but for tlie merger model. 
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is based on the lines eon tained in the CD23 eompilation 
of iKurucz fc Bel] ()1995l) . To account for higher ioniza- 
tion at early times, we added the ionization stages VI and 
VII for Sc to Ni, leading to a total of ^5 x 10'^ atomic 
lines. 

Both our models yield a B-band peak magnitude of 
— 19.0, roughly in agreement with that observed for 
SN 201 Ife. Their rise times, however, differ: while NlOO 
reaches i?-band maximum after 16.6 d, the merger takes 
20.8 d (further parameters of our synthetic light curves 
are given in Table [T]). Thus, neither of the models gives 
a perfect match to the light curves of SN 2011fe but both 
are sufficiently close to warrant further investigation. 

In Fig. [3] we compare synthetic spectra from our mod- 
els with flux-calibrated spectra of SN 201 Ife taken by 
the SNfactory collabo ration with the SNIFS instrument 
(jAldering et al.l[2002[ ) on the University of Hawaii 2.2 m 
telescope on Mauna Kea. To our knowledge this is the 
first direct comparison of consistent three-dimensional 
SN la models with a spectrophotometric time-series. 
Overall, the spectra of both explosion scenarios repro- 
duce the main features of the observed spectra and the 
flux level reasonably well (note that these are not fits but 
predictions from "first-principle" models and that abso- 
lute fluxes are compared). In detail, however, there are 
problems in both models. 

The predicted absorption features are blue-shifted with 
respect to the observations. Although this effect is 
stronger in NlOO, it is also visible for the merger, in- 
dicating too high ejecta velocities in the models. Since 
the mass of the exploding object is very different in the 
two cases, the nuclear energy release is likely too high in 
both explosion processes. A potential way to cure this 
problem is to increase the oxygen abundance in the pro- 
genitor WDs at the expense of carbon thus increasing 
the average nuclear binding energy of the fuel. 

While NlOO is only marginally too bright at the early 
epochs, the merger is clearly too faint. This corresponds 
to the shorter/longer rise time of the respective model 
(see Tab. [T|) compared to SN 201 Ife. Around B-hand 
maximum at ^18 d the merger compares favorably to the 
observed spectra. The flux level and the overall shape of 
its synthetic spectra match the data better than those 
of NlOO. After maximum the agreement with the obser- 
vations deteriorates for both models, although the effect 
is more drastic for NlOO. In particular, the models fail 
to reproduce the spectral features between 5000 A and 
6000 A. Moreover, both models become redder faster 
than the observation. Again, this trend is more pro- 
nounced in NlOO, but is also visible in the merger. 

Figs. [1] and [2] show that iron-group elements (IGEs) 
extend signiflcantly beyond velocities of 10, OOOkms"^ 
in NlOO but also in some direction s in the merger. The 
W7 model of iNomoto et all (|1984f ). which is known to 



reproduce observables of SNe la well, does not contain 
IGEs at such high velocities. However, they are reported 
in abundance tom ographies of the normal SNe 20 02bo 
(iStehle et al.l[2005h a nd 2004eo (jMazzali et al.l[2008l ) and 
iNugent et al.l ()2011b( ) identify iron in the earliest spectra 
of SN 2011fe at velocities of 16,000kms'^ Our syn- 
thetic spectra do not show mismatches with observed 
lines that can be directly attributed to high-velocity 
IGEs. It is possible, however, that they contribute to 



the fast reddening of the models. As for the high ejecta 
velocities, a decreased carbon/oxygen ratio in the ex- 
ploding WDs may alleviate this problem. 

An important difference between the two models 
is also visible from Figs. [1] and [2l NlOO is — apart 
from small-scale ani sotropics — roughly spherical (see also 
iBlondin et al.|[20Tl|) . In contrast, the merger shows pro- 
nounced large-scale asymmetries. This is reflected in a 
strong viewing-angle dependence in its spectra which in- 
creases after maximum due to the growing asymmetry of 
the inner ejecta for smaller radii. As demonstrated in the 
lower plot of Fig. [3] individual line-of-sight spectra from 
the merger reproduce the observation at least as well as 
the angle-average. Nevertheless, the high level of asym- 
metry in the merger could be in conflict with the observed 
spectral homogeneity of normal SNe l a and the low level 
of co ntinuum polarization observed (jWang fc Wheelen 
|2008| ). Currently our RT simulat ions do not include po- 
larization. Note, however, that iSmith et al.l (|2011[ ) re- 
port a continuum polarization of 0.2-0.4% for SN 2011fe 
which they interpret as a sign of persistent asymmetry 
in the last-scattering surface. 

The fact that the merger reproduces the observed spec- 
tra better than NlOO at maximum light (and later) sug- 
gests that the chemical structure of its deeper ejecta, 
which dominate the spectrum formation at this epoch, 
is closer to that of SN 201 Ife. Since neither of the mod- 
els matches the optical data of early epochs perfectly, 
our comparison gives slight preference to a WD merger 
scenario over a delayed detonation in a Chandrasekhar- 
mass WD as an explanation for this object. But as 
there are major shortcomings in both models, a defini- 
tive conclusion cannot be drawn. Observables other than 
maximum-light spectra may, however, have more dis- 
criminating power. We discuss promising possibilities in 
the following. 

4. LATE-TIME OBSERVABLES 

While optical data taken before and around peak 
brightness probe predominantly the outermost layers, 
observations at later epochs are sensitive to the core of 
the ejecta. Since this is the region where the differences 
between the two models considered here are most pro- 
nounced, late-time observables are a very useful diagnos- 
tic tool. 

A fundamental difference between our two explosion 
scenarios is the density at which the material is burned 
in the thermonuclear combustion. Due to the high cen- 
tral density of the Chandrasekhar-mass WD, substantial 
burning proceeds on thermodynamic trajectories with 
peak densities above 2 x 10* g cm~'^ in NlOO (especially in 
its deflagration phase), whereas in the violent merger all 
the burning occurs at peak densities below 2 x 10* g cm~^. 
This leads to vastly different degrees of neutronization 
in the ashes due to electron capture reactions, result- 
ing in higher abundances of stable IGEs (in particular 
54pg^ 56pg^ ^^^ ssj^j) foj. |>^jQQ ^jjjpjj should be reflected 

in the presence of Ni lines in spectra from the nebular 
phase (jMaeda et al.ll2010bt ICerardv et al.ll2007| ). More- 
over, our models differ significantly in composition and 
geometry of the innermost ejecta. While the ^^Ni distri- 
bution is roughly spherical for NlOO, most of the ^^Ni is 
off-center in the merger. Asymmetric and shifted lines 
can be expected here and could correspond to the effects 
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Figure 3. Spectral evolution of our delayed detonation NlOO (left) and merger model (right) from 6 to 27 days after the explosion. The 
angle-averaged spectrum is plotted in black while 25 spectra for representative viewing angles are shown in gray (the variability with viewing 
angle of the earliest spectra is dominated by Monte Carlo noise in both models). For comparison the observed spectra of SN 2011fe are 
over-plotted in red assuming an explosion date at Augu st 23.7 (MJD 55796. 7; [Nugent et al. 2011b). The observations were corrected for 
Galactic reddening assuming E{B — V)Gal = 0.009 mag l lSchlegel et al.iri998l ) and dc-redshifted according to a heliocentric radial velocity 
iijiel = 241kms~^ given bv lde Vaucouleurs et all 1 119911) . Reddening from the host is negligible. The bottom panel compares the observed 
spectrum of SN 201 Ife near B-band maximum (red) with synthetic spectra from the merger corresponding to three different viewing angles 
(blue colors) and the angle-average (black). 
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discussed by iMaeda et alj (l2010bl ia[). In the merger, the 
detonation of the secondary WD at low densities pro- 
duces copious amounts of oxygen in the innermost cjecta. 
Potentially, this could lead to strong [01 ] AA6300,6364 
emission at late times (|Kozma et al.l 120051 ) . which is not 
observed in SNc la. The efficiency of [01] emission, how- 
ever, depends strongly on the contamination with other 
elements and the ionization state of the oxygen- rich zone. 
For strong ionization (as might be the case in our low- 
density core) , the presence of oxygen in the core would 
pose no problem for the merger. Clearly, nebular spec- 
tra contain important information, but a firm conclusion 
can only be drawn upon detailed three-dimensional mod- 
eling of the late-time RT, which is beyond the scope of 
this letter. 

There is, however, a more direct and perhaps more 
easily studied effect. Due to the low central densities of 
the two sub-Chandrasekhar mass WDs. all of the IGEs 
in the merger are synthesized in either a-rich freeze-out 
from nuclear statistical equilibrium or in incomplete Si- 
burning. In contrast, much of it is produced under "nor- 
mal" freeze-out conditions in the defiagration phase of 
NlOO (Thielemann et al. 19861 place the dividing line be- 
tween a-rich and "normal" freeze-out at ~2 x 10^ gcm~'^ 
for explosive burning of C -1-0 material in SNe la). This 
leads to a higher abundance of ^^Co — an isotope mainly 
synthesized in t he "normal" freeze-out a nd in incomplete 
Si-burning (e.g. lThielemann et al.lll986[ ) — in the ejecta of 
the Chandrasekhar-mass WD explosion than in those of 
the merger. 

Such different isotopic ratios in the IGEs affect 
the shapes of the predicted late-time light curves 
(jSeitenzahl et a l. 2009b). Starting at --800 d after the ex- 
plosion, the leptonic light curves that assume full trans- 
parency to 7-rays and pure leptonic heating of the ejecta 
will be increasingly powered by the decay of isotopes 
other than ^^Co. This is illustrated in Fig.Sl At -1000 d 
after the explosion, the decay of ^^Co to ^^Fe, which (in 
—80% of all decays) emits internal conversion electrons, 
starts to dominate the light curves. Later, the decay of 
^^Fe, which is mainly synthesized as ^^Co, to ^^Mn (a 
ground-state to ground-state transition followed by the 
emission of Auger electrons) contributes significantly and 
eventually dominates the radioactive energy generation. 

The leptonic light curve of NlOO will fall off more slowly 
than that of the merger. For example, the decrease 
in combined leptonic energy production from 900 d to 
1400 d (1900 d) corresponds to a dimming by 2.25 (3.20) 
magnitudes for NlOO and by 2.65 (3.87) magnitudes for 
the merger (see Table [1]). 

Thus, a measurement of the late light curve de- 
cline rate would distinguish between an explosion of 
a Chandrasekhar-mass WD (which in any scenario re- 
quires some pre-expansion in a deflagration stage) and 
alternative models based on detonations in low density 
material — such as mergers of WDs. For this, neither the 
correct distance to the object nor the exact '^^Ni produc- 
tion have to be known. Note, however, that the light 
curves shown are idealized cases assuming that the lep- 
tonic energy production rates can be directly translated 
into UVOIR light curves (which may be precluded by 
effects such as the infrared catastrophe, "frozen-in ion- 
ization", CSM interaction, leptonic losses, etc.). 
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Figure 4. Leptonic luminosity as a function of time after the ex- 
plosion for NlOO (black bold line) and the merger model (red bold 
line). The dashed, dash-dotted and dotted lines give the contribu- 
tion to the leptonic luminosity due to the decay of selected isotopes. 
Thin solid lines give the luminosities including decay X-rays. 

5. CONCLUSIONS 

Although the nearby SN 2011fc offers a unique oppor- 
tunity to scrutinize explosion models, at present a clear 
preference of one scenario over the other is hard to es- 
tablish. We therefore discuss two models that are very 
distinct in the explosion characteristics and in the re- 
sulting structure of the ejecta — a delayed detonation of 
a Chandrasekhar-mass WD and a merger of two WDs 
with a total of 2 M©. 

Comparing with early and near-maximum optical spec- 
tra, both scenarios reproduce the main features but the 
merger is slightly preferred because it provides a better 
match to the observations around peak brightness. There 
are, however, shortcomings in other aspects — such as the 
too long rise time — and therefore the working hypothesis 
of SN 2011fe resulting from a merger of two WDs requires 
additional confirmation. 

As shown here, alternatives to early-phase optical data 
may have more decisive power. At very late epochs nu- 
cleosynthetic effects lead to different characteristics in 
the photometric evolution that may allow us to discrim- 
inate between explosion models. This, however, requires 
true bolometric measurements, and it is unclear at which 
wavelengths the maximum emission occurs at those late 
epochs. Observations of SN 201 Ife will help to clarify 
this issue and thus multi-wavelength monitoring of this 
object is essential. If the maximum emission falls into the 
optical range, a clear distinction between explosion mod- 
els (that can then be related to progenitor scenarios) will 
be possible from photometric measurements at >1000d. 
Thanks to its proximity, these observations should be 
feasible for SN 2011fe. 
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